Exciton Hierarchies in Gapped Carbon Nanotubes 
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We present evidence that the strong electron-electron interactions in gapped carbon nanotubes 
lead to finite hierarchies of excitons within a given nanotube subband. We study these hierarchies 
by employing a field theoretic reduction of the gapped carbon nanotube permitting electron-electron 
interactions to be treated exactly. We analyze this reduction by employing a Wilsonian-like numer- 
ical renormalization group. We are so able to determine the gap ratios of the one-photon excitons 
as a function of the effective strength of interactions. We also determine within the same subband 
the gaps of the two-photon excitons, the single particle gaps, as well as a subset of the dark exci- 
tons. The strong electron-electron interactions in addition lead to strongly renormalized dispersion 
relations where the consequences of spin-charge separation can be readily observed. 
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Gapped carbon nanotubes are the subject of intense 
experimental [l|-|6[ and theoretical interest [7Hl7| due to 
their possible application as opto-electronic devices [l| 
as well as providing a particularly clean realization of 
str ong ly correlated electron physics in low dimensions 
[9|,ll7|. Both interests converge in the study of the tubes' 
excitonic spectra. This spectra, due to the effects of a 
strong Coulomb interaction in one dimension, is strongly 
renormalized from what would be expected from the un- 
derlying band structure of the tube. Its computation 
therefore requires a non-perturbative approach. 

The favoured theoretical approach to studying the ex- 
citonic spectra of carbon nanotubes is the use of a Bethe- 
Salpeter equation combined with first principle input to 
approximate the particle and hole wavefunctions that 
form the excitons [lOJ, luJ, [iSl, ll5 1 . This approach has been 



particularly valuable in that it has allowed a quantitative 
description of aspects of excitonic physics, in particular 
the magnitude of the excitonic gap for the lowest lying 
excitons in a given subband of the nanotube. In this 
letter we step away from trying to describe quantitative 
details of the excitons and instead focus on more qualita- 
tive features. We do so using an approach that combines 
a field theoretic reduction of the nanotubes identical to 
that used to study Luttinger liquid behaviour in metallic 
nanotubes with a numerical renormalization group that 
enables one to study the effects of gapping out a multi- 
component Luttinger liquid. In the process we find in a 
regime of strong electron-electron interactions a number 
of new qualitative features to the excitonic spectra (see 
Fig. 1). 

First and foremost we find that a given subband of 
the nanotube has a finite number of optically active one- 
photon excitons, A-'^"'*, where the multiplicity depends 
on the strength of the tube's screened Coulomb interac- 
tion. While a simple picture of excitons as analogs of the 
excited states of a hydrogen atom [7| yields an infinite hi- 
erarchy of excitons, we argue this series is truncated to at 
most three. In order to understand the excitons' binding 



energy, we also study the single particle spectrum of the 
tube. The single particle gap is strongly renormalized 
by interactions and can be many multiples of the bare 
bandgap. Consequently the lowest lying two-excitation 
continuum, depending on strength of the Coulomb inter- 
action, may be a particle-hole continuum or it may be a 
two-exciton continuum. This approach treats all excita- 
tions of the nanotube on the same footing. It thus con- 
flates the difference between an exciton and a bi-exciton. 
Indeed interactions can be sufficiently strong that exci- 
tons should then be thought of equally as particle-hole 
bound states or bound states of two other excitons. 
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FIG. 1: A schematic indicating the position in energy of the 
various excitons together with two-excitation continua arising 
from our treatment at larger values of the effective interaction 
(A'^^ > 5). The bracketed numbers are estimates computed 
at Kc+ = 1/7 of the relative absorption strength (relative to 
the strength of A^"'^) of the three one-photon optically active 
excitons, A^""-^-^. 



Our results are applicable both to semiconducting nan- 
otubes and metallic nanotubes which have been turned 
semiconducting by the application of an axial magnetic 
field l8|. In the latter case where the gap is tunable and 
can be made much smaller than the bandwidth, we ex- 
pect our approach to yield results that are quantitatively 
accurate. For tubes that are semiconducting due to band 
structure, our results within a one subband approxima- 
tion are most reliable for large radius tubes although 
we believe the qualitative features that we find to apply 
equally well to smaller radius tubes. 

We focus on a single subband of a carbon nanotube. 



To describe this subband at low energies we introduce 
four sets (two for the spin, a, degeneracy and two for 
the valley, a = K,K', degeneracy) of right (r = +) 
and left (r = —) moving fermions, ipraa- The Haniil- 
tonian governing these fermions can be written as H — 
J dx{Hkin + Hgap) + Hcouiornb ■ Hki„ + Hgap together give 
the non- interacting band dispersion, e^(p) = Vqp'^ + Aq: 

Hkin = ~iV0'4'laaC>x-ipracr] Hgap = Aalpl^^lp^racr, (1) 

where vq is the bare velocity of the fermions and re- 
peated indices are summed. For the Coulombic part of 
the Hamiltonian we include only the strongest part of 
forward scattering: 



Hco 



dxdx' p{x)Vo{x — x')p{x'), 



where p{x) = J2raa'^taaix)'^racrix). Wc also ignore 
backscattering interactions arising from Coulomb inter- 
actions. In the presence of a gap term, these interactions, 
as they are marginal in the RG sense, are less important. 
But as they lift the model's SU(4) symmetry, they will 
split the lowest lying 15-fold multiplct of optically dark 
excitons. 

To study the full Hamiltonian, wc treat Hgap as 
a perturbing term (albeit one to be treated non- 
perturbativcly) of Hq = J dx{Hkin) + Hcouiomb- The 
latter terms are nothing more than the Hamiltonian of 
a metallic carbon nanotube. In excitonic language, we 
thus think of Hgap as a confining interaction on top of 
the metallic tube. The advantage of doing so is that we 
are able to treat Coulomb interactions analytically ex- 
actly at the beginning of the computation. We do so 
using bosonization. 

If we bosonize Hq in terms of chiral bosons (jjraa by 

wri ting ipraa- ~ exp(i0rao-), we arrive at a simple result 

17|,|l8|. The theory is equivalent to four Luttingcr liquids 

described by the four bosons 9i, i = c±,s± (and their 

duals (j)i) 



Ho= fdxY,^ (M^xC^^)^ + Kr\dM' 



(2) 



The four bosons diagonalizing Hq are linear combina- 
tions of the original four bosons and represent an effec- 
tive charge-flavour separation where 6c+ = X^rao- ''"'Praa 
is the charge boson and the remaining three bosons reflect 
the spin, valley, and parity symmetries in the problem. 
The charge boson is the only boson to see the effects 
of the Coulomb interaction. Both the charge Luttinger 
parameter Kc+ and the charge velocity Vc+ ~ vq/Kc+ 
are strongly renormalized. In particular for long range 
Coulomb interactions, Kc+ takes the form 



A'c+ = I 1 + 



8e^ 



L 



ttkHvo 



(log(:^) + co) 



'2TrR 



-1/2 



(3) 



where k is the dielectric constant of the substrate, L is 
the length of the nanotube, R is the tube's radius, and 
cg is an 0(1) numerical constant dependent on details of 
the wrapping vector (n, m) [16lll7l|. In typical nanotubes 
Kc+ can take on values in the range of ^ .2. With an elec- 
tric field applied transversely to the tube, Kc+ can easily 
be made as small as .1 [l6|. The remaining Luttinger 
parameters, Ki, i = c-,s± retain their non-interacting 
values, 1, and so their velocities, Vi = vq go unrenor- 
malized. Including the full forward scattering leads to a 
small (upwards) renormalization (on the order of 10%) 
of these parameters which we will not consider here [17[ . 
These rcnormalizations, however, will not change any of 
our results qualitatively. 



Under bosonization, Hgap becomes 



Hgap — 



4Ao 



(J|cos(|) + nsin(^)), 



(4) 



where Ao = Ao{A/vc+Y^~^'--+^^^ and A is the bandwidth 
of the tube. This gap term is highly relevant. Rewriting 
the gap term in this fashion already gives us important 
generic features of the excitonic spectrum. Firstly, as a 
perturbation of Hq, Hgap has the anomalous dimension, 
3/4 -I- Kc+/4:. In turn this implies the full gaps, {Aq}q,, 
satisfy the scaling relation 



A. = A(^)^/(^-^^.)/„(if,^. 



(5) 



where {fa} a is a set of a priori unknown dimensionless 
functions of Kc+. We will see that for the excitons, Jexc 
is a relatively weak function of K^^ while for the sin- 
gle particle excitations, fsp depends strongly on A'c+. 
This scaling relation tells us immediately how exciton 
gaps scale between subbands. If we take a large radius 
tube where the gap of the n-th subband gaps scale as 
Ao(n, i?) ^ n/R, the expected gap ratio between the 
first and second subband, is 2^/^^"-^'^+) ^ 1.78 (this is in 
rough correspondence to that reported in ^ and so pro- 
vides a straightforward explanation of the ratio problem 
12|). Secondly, Hgap takes the form of a perturbation of 
a generalized sine-Gordon model (involving four bosons 
instead of one). We thus expect that excitons (within a 
given subband) will come in hierarchies whose size is de- 
termined by the value of Kc+ very much like the number 
of bound states in a sine-Gordon model, a perturbation 
on a free boson of the form cos(/30), is determined by the 
parameter (3. 

Having completely characterized Hq , we study the con- 
fining effects of Hgap using the truncated conformal spec- 
trum approach (TCSA) [l9| combined with a Wilso- 
nian renormalization group modeled after the numeri- 
cal renormalization group (NRG) used to study quan- 
tum impurity problems [20[. This methodology permits 
the study of arbitrary continuum one dimensional theo- 
ries provided one can write the theory as the sum of a 



gapless theory, in this case Hq, plus a perturbation, here 
Hgap- It is flexible enough that it can play the same role 
as DMRG does for lattice models (although as it is an 
NRG, it accesses matrix elements more easily). Among 
the first uses of the original TCSA, unequipped with an 
NRG, was to study a quantum Ising model in a magnetic 
field [l9[ . The model we study here, in terms of compu- 
tational complexity, is equivalent to eight coupled Ising 
models. Because of this added complexity, this model 
could not be studied using the TCSA without the ac- 
companying renormalization group. Additional details 
of this approach as applied to the problem at hand may 
be found in Section I of Ref. 21 1. 



Our results for the gaps as a function of Kc+ are pre- 
sented in Fig. 2. Details of the analysis of the numerical 



data have also been relegated to Ref. |2l| . We emphasize 
again that these excitations are all found within a single 
subband of the tube. We first consider the one-photon 



optically active excitons (labelled A 



here u refers to the antisymmetry under tt— rotations 
about the centre of a carbon hexagon - in bosonic lan- 
guage antisymmetry under this rotation corresponds to 
antisymmetry under 9c+.s+ -^ ^&c+.s+)- We see that as 
Kc+ decreases (i.e. the effective Coulomb interactions be- 
come stronger) the number of one photon excitons grows. 
At Kc+ > 1/5 the number of such excitons goes from one 
to two, while at Kc+ > 1/6 the number of such excitons 
goes from two to three. But thereafter this number sat- 
urates. 

There are two limits on stable excitons in the sys- 
tem - the beginnings of both the particle-hole contin- 
uum (marked in dark green in Fig. 2) and a two-exciton 
continuum (marked in maroon in Fig. 2). Necessarily ex- 
citon gaps cannot cross over these boundaries lest a decay 
channel opens to the exciton. The bottom of the particle- 
hole continuum is twice the single particle gap, 2A'*p. As 
A^P is a strong function of Kc+ (going as A**^ ~ K^_^ ) , at 
small Kc+ there is room for additional excitons. This is 
very much like the sine-Gordon model where as (3 goes to 
zero, the number of bound states goes to infinity. How- 
ever for an exciton to be stable it must also fall below 
a two-exciton continuum. This two-exciton continuum, 
formed from the lowest single-photon exciton, A^"'^, and 
the lowest two-photon exciton, A^^'^, is relatively insen- 
sitive to Kc+ and serves to provide a bound on the total 
number of possible excitons at small Kc+- 

In bosonized models considerable intuition can be had 
by considering the classical limit and asking what are 
the classical field solutions that correspond to the exci- 
tons/bound states. This is true here as well. The first 
one-photon exciton, A^"'*, corresponds to a breather- like 
solution of the field equations governing Hq + Hgap- This 
solution corresponds to 9c+ interpolating between to 
Att and back to while the other three bosons remain 
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FIG. 2: The exciton and single particle energies as a function 
oiKc+. 



1, 2, 3 where fixed at 0. It takes the form 



7c+ 



{x,t) = 4 tan 



\/l — Ld^ cos{ujt) 
Lo cosh(\/l — oj'^x) 



(6) 



where a; is related to the (classical) energy of the 
breather. We can estimate the energy corresponding 
to quantizing this classical solution via mean field the- 
ory. Replacing Hgap by 4Ao/7rS^ cos(0c-(-/2), where S = 
(cos(0//2)) with Of — 6c-,ds+, or 0s_, we reduce the 
model to a pure sine-Gordon theory. A^"'^ is then equal 
to the gap of this theory's first bound state [23[ 



lu,l 
^MFT 



= 2(2Ao^ 



r(i 



^'c+ 



)^^ 2^(|)sin(f) 



(7) 



where a ~ \/AK'^j^ and ^ = ct/(2 — ct). S is 
determined self-consistently by replacing Hgap 
with 4Ao7r"^S2(cos(6'c+/2))cos(6'//2) and using 
22| to express S in terms of the effective cou- 
pling constant 4Ao7r~^S^' 



quently find that as K^ 



c+ 



A 



MFT 



/(A(Ao/A)4/(5 



-K, 



cos(6'c+/2)). We conse- 
varies between 2 and 10, 
varies from 3.1 to 4.1 in 



reasonable agreement (if shifted upwards) with the 
non-perturbative NRG. 

The second one-photon exciton, A^"'^, has a gap a 
little more than twice that of A-'^"'^ and corresponds to 
roughly the putative location of A^"'^, the first exciton 
in the second subband (which we have argued, by scal- 
ing, is found at ^ 1.8A^"'^). If we were to take into 
account intersubband Coulomb interactions, hitherto ig- 
nored, A^"'^ and A^"'^ would mix leading to a set of 
hybridized states. It is thus conceivable that reports of 
observations of the exciton A^"'^ are in fact seeing A^"'^ 
or at least some hybridized combination of the two. 

Beyond the one-photon excitons. Fig. 2 also displays 
the gaps of a number of other excitons. It shows the 
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FIG. 3: The dispersion of the first two excitons A^"'^ and 
A^"'^ for Kc+ = 1/5. Dashed lines in the small, but finite, 
momentum regions connecting data points are guides to the 
eye. 



single existing two-photon exciton, A^^'-'^ (here g cor- 
responds to excitons which are symmetrical under tt- 
rotations about a centre of a carbon hexagon). Unlike 
the one-photon excitons, there is only a single stable two- 
photon exciton regardless of the value of Kc+- Lying 
just below A^"'^ are found a set of 15 dark excitons with 
energy A'^"'"'^. These excitons subsume the dark triplet 
excitons but are of greater degeneracy because they not 
only carry spin quantum numbers but valley quantum 
numbers as well (dark excitons with non-trivial valley 
quantum number were termed i^-momentum excitons in 
Ref. @). As the full symmetry of H is SU(4), the dark 
excitons transform as SU(4)'s adjoint. At twice the en- 
ergy of A'^'^^'' begins a two exciton continuum into which 
any two-photon exciton will decay. As A'^^'^ is just be- 
low 2A'^°'''^, we understand this two exciton continuum 
essentially precludes any two-photon excitons with en- 
ergy greater than A^^'^. 

By examining the momentum dependence of the exci- 
tons' dispersions, we can see the consequences of having 
two velocities, charge, Vc+, and flavour vq — Vc-,Vs±, 
in the system. The above semiclassics suggests that the 
initial dispersion of A^"'^, £^^"'^(p), will behave approx- 
imately as ((A^"'^)^ + p^v-^^y/^ (a consequence of it in- 
volving 9c+ degrees of freedom alone). While we can- 
not write down the classical solution corresponding to 
A^"'^, it lies primarily in the spin/flavor sector and so 
will disperse governed by the much smaller velocity, vq. 
While A^"'^ > A^"'^, the much larger charge velocity 
means that the dispersion of A-'^"'^ will intersect A^"'^ 
at a relatively small momentum, pint = (((A^"'^)^ — 
(A^"'^)2)/(^2^-'(;o))^/2. This intersection means the two 
states will hybridize. As a consequence of this hybridiza- 
tion, for momenta larger than pi„t, A^"'^(p > pint) will 
disperse with an effective mass determined by vq whereas 
A^'^''^{p > Pint) will now vary according to a velocity 
much closer to the larger Vc+. 

We observe this general behaviour in our numerical 



analysis. Plotted in Fig. 3 are the dispersions of these 
two excitons for K~^ = 5. We see that E^"-'^{p) in- 
creases extremely rapidly for momenta up to p ^ A/2uo 
but thereafter levels out, even decreasing. This increase 
is greater than would be predicted from the relativistic 
dispersion, {{A^^-^y + p'^v^_^_y^^ , suggesting that charge 
and spin, even at the smallest momenta, are intricately 
linked, leading to an even stronger renormalization of the 
effective mass. If however we combine the renormaliza- 
tion due to Vc+ with the logarithmic correction p'^log(p), 
to the exciton self-energy predicted in Refs. [12|, |l3[, we 
can understand the observed behaviour. Complementar- 
ily, E^^''^[p) initially displays a weak momentum depen- 
dence but then increases with a velocity approximately 

Vc+ll for p > Pint- 

Finally we touch upon the excitonic signal in absorp- 
tion spectra. The strength of this signal is proportional 
to the imaginary part of the current-current correlator. 
It is straightforward to compute the necessary matrix el- 
ements of the current operator within the framework of 
TCSA+NRG. We find, as indicated in Fig. 1, that A^^'^ 
has by far the strongest signal. While A^"'^ should be 
visible in an absorption spectra, A^"'"^, with a weighting 
1/1000 of A^"^! is effectively dark although could be de- 
tected indirectly through its effects on the temperature 
dependence of the excitonic radiative lifetime [l3| . 

In conclusion, we have argued for, using a fully many- 
body, non-perturbative approach, new features in the 
excitonic spectrum of gapped carbon nanotubes. We 
have shown that with a given subband, optically ac- 
tive excitons appear in finite hierarchies. We have also 
shown that the single particle sector experiences ex- 
tremely strong renormalization so swapping the begin- 
ning of the particle-hole continuum with the two-exciton 
continuum. Further work will include studying the ef- 
fects of backscattering interactions in the tubes (not 
treated here) together with intersubband interactions. 
More generally, we have demonstrated a robust method 
able to treat continuum representations of arbitrary one- 
dimensional systems. 
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SUPPLEMENTARY MATERIAL: EXCITON 

HIERARCHIES IN GAPPED CARBON 

NANOTUBES 

Description of NRG treatment of TSA as applied to 
gapped carbon nanotubes 

We study the confining effects of Hgap using tlie trun- 
cated conformal spectrum approacli. As orginally de- 
veloped [19[, the TCSA has been used extensively (for 
a sample consider 
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-30[). However here we combine 
it with a Wilsonian renormalization group modeled af- 
ter the numerical renormalization group (NRG) used to 
study quantum impurity problems [20|. With quantum 
impurity problems, the NRG exploits the idea that sites 
on the Kondo lattice far away from the quantum im- 
purity are less important than sites close to the impu- 
rity in determining the physics. We have a similar no- 
tion here as the perturbation Hgap is relevant in the 
RG sense (the anomalous dimension of the operator, 
3/4 + A'c+/4, is less than 2) and so the eigenstates of 
the unperturbed Hamiltonian, Hq, lowest in energy are 
most important in determining the physics of the full 
model. We are thus justified in initially truncating the 
spectrum oi Hq = J dx{Hkin)+Hcouiomb to a finite num- 
ber of states (for the computations chosen here, we ini- 
tially truncate at 2500). To be able to do so, we study 
the theory at finite size, Rnbg- In making the system 
finite, the spectrum becomes discrete with level spacing 
2t: /Rnrg. However even at finite size, the conformal 
nature of Hq means we have complete control over the 
spectrum (which can be thought of as a four-fold tensor 
product of states coming from Hilbert spaces of compact- 
ified bosons). 

The Hilbert space of the theory can be described more 
properly as follows. For a single boson, with compacti- 
fication radius, R, the Hilbert space is built from a set 
of highest weight states denoted by \n,m) where n gives 
the center of mass momentum, n/ R, of the state, and m 



denotes the winding number (or charge) of the state [31 
\n,m) is related to the true vacuum state, |0), by 



\n,m) 



e'RHo)+^Ho)\Q^^ 



where (p is the boson dual to 9. Descendant states may 
be created by acting on the |n, to) with a set of oscillator 
modes a_„, a_„, n > 0, that come from mode expanding 
the boson. For our multibosonic system with bosons, 9c±, 
0s±, the set of highest weight states are formed according 
to 



'E, 



-e.(o)+»E, j^0j(o) 



y rii mod 4 = 0; >, "^i mod 4 = 0; 

i i 

y niirii mod 4 = 0; 



{rii + rii^i) mod 2 = 0; {rrii + 77ii+i) mod 2 = 0. (8) 

These particular conditions upon the highest weight 
states may be derived by considering the original bosons 
that arise in bosonizing the subband fermions (i.e. be- 
fore the transformation to the charge and spin bosons). 
Descendant states are again formed by acting on these 
states with the modes ai,_„, ai._„, n > 0, of the bosons. 
We note that we can show this Hilbert space leads to a 
modular invariant partition function indicating we have 
a well defined theory [31| . The action of Hq upon a state 
of the form 



\E') 



Ho\E') = 
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-1/12) 



-Y,v,iY,k,i)W^), 



1=1 



(9) 



where the Kj are the four Luttinger parameters. 

Having described the Hilbert space, we now label the 
discretized spectrum, \E^), ji^j), • • • where the spectrum 
has been ordered in terms of increasing energy. We then 
truncate the spectrum keeping the first A^ + A states of 
the unperturbed theory, |£'i), • • • , I^^at+a) (where N = 
1500 and A = 1000 were typically chosen for this study). 
We further exploit our control over the unperturbed the- 
ory by computing all the matrix elements of Hgap, i.e. 
{E^\Hgap\Ej). This knowledge then allows us to repre- 
sent the full Hamiltonian as a (A^-|~ A) x {N + A) matrix 
which we then diagonalize. If Hq were a simple theory, 
this procedure would be enough to obtain highly accu- 
rate results. In studying the critical transverse field Ising 
model perturbed by a parallel magnetic field, Ref. llOlwas 
able to keep only 39 states and yet still obtain the spec- 
trum of the theory to an accuracy of a few percent. How- 
ever a single Luttinger liquid can roughly be thought of 
(in terms of complexity of the Hilbert space) as a tensor 
product of two Ising models. Thus Hkin + Hcouiomb, be- 
ing composed of four bosons, can be thought of as equiv- 
alent in complexity to a tensor product of eight Ising 
models. This eight fold complexity means that a sim- 
ple truncation can no longer hope to describe correctly 
the physics. This is where the numerical renormalization 
group (NRG) is crucial, which we now briefly describe. 

The NRG begins by taking the states arising from 
the diagonalization of the {N + A) x {N + A) ma- 
trix and ordering them in terms of increasing en- 
ergy, |£^i), • • • , I-EjV+a)- We then discard the top A 
states. We re-form the truncated Hilbert space by 
taking the lowest N energy states remaining from the 
diagonalization and add to them the next A states 



from the unperturbed spectrum. We thus obtain a 
truncated Hilbert space of the same size (A'^ -|- A): 
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We re-form the 



\El),- 

Hamiltonian in this new basis, rediagonalize, and then 
repeat the procedure until convergence. In this way, we 
allow the high energy portion of the unperturbed spec- 
trum to influence the low energy sector in bite size, nu- 
merically manageable (A) chunks. 

In our study of the exciton gaps, we typically iterate 
the NRG up to 70 times (and so with a A = 1000, we take 
into account 70000 states of the unperturbed theory). 
While the results are not completely converged, we allow 
an analytic RG to take over from the numerical RG. At 
sufficiently high energies, the flow of a gap under the RG 
is governed by a simple one loop equation [20[. We can 
thus flt the flow equation to our numerical data at high 
truncation energies and extrapolate to infinite truncation 
energy. For single particle gaps however, we found we 
needed to iterate the NRG many more times than for the 
excitons. Because of the large charge velocity (at least at 
small values of Kc+ ) , we needed to take into account on 
the order of 1200000 states from the unperturbed theory. 
If we were to use the NRG as described in Ref. [20|, 
this would mean iterating the NRG over 1000 times, a 
too great of numerical burden (the computational time 
scales as the number of iterations to the third power). 
We have thus adapted the NRG in this case. The great 
majority of the 1200000 states that the NRG needs to 
run through in the single particle case do not influence 
the end result. We thus run a 'preselection' of the most 
important states and run the NRG on those. This works 
straightforwardly as follows. 

To run the preselection, we run the NRG as initially 
described for M iterations (and so account for the first 
N+MA states of the unperturbed theory). We have then 
obtained a reasonable approximate (though not yet fully 
converged) of the lowest lying eigenstates of the theory. 
Let us for the sake of specificity focus on the ground state, 
I^Gs)' obtained after these M-iterations. Rather than 
blindly run the NRG further, we use these eigenstates 
to provide a measure of which of the yet to be iterated 
through states are important. Suppose we want to iterate 
through a further KA states. We first compute their 
overlap with {Eq^) in second order perturbation theory, 
i.e. we compute 
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We take the Ti 's as a measure of how strongly the state 
\E%^j,j^j.^) will mix into \E^^) under the NRG. Using 
this measure, we then throw away states, l-E-jv+MA+i) 
with Ti below some cutoff, A. We then take the remain- 
ing states that have survived this cutoff procedure and 
run the NRG on those as before. 

Typically in applying this procedure to the single par- 
ticle gaps, we have used M = 70, A = 1000, K up to 
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FIG. 4: The value of A^"'^ for Kc+ = 1/5 as function of the 
system size, Rnrg, and the cutoff energy, Ec- 



1100, and A = 2 x 10"'*^. After applying the trunca- 
tion to the KA states, we were typically left with 70000 
states. Running through these 70000 states is then done 
in several hours rather than the several weeks it would 
take to run through the original KA states. 

In analyzing the Hamiltonian we are careful to mini- 
mize the necessary computational burden by respecting 
the symmetries of the full Hamiltonian. The full Hamil- 
tonian possesses a C/(l) x 5'[/(4) symmetry. The C/(l) 
arises from the (valley symmetric) charge degrees of free- 
dom while the SU{A) governs the flavour degrees of free- 
dom (composed of valley antisymmetric charge degrees of 
freedom plus spin degrees of freedom) . The Hamiltonian 
also possesses Z2 degrees of freedom which arise by trans- 
forming 2 of the 4 bosons, 9i, via 9i — > —6i. The excitons 
that are optically active are all S'C/(4) scalars and so we 
run our numerics in this portion of the Hilbert space of 
the full Hamiltonian. Similar the single particle excita- 
tions transform as a four dimensional representation (the 
fundamental under SU(4)) while the one branch of dark 
excitons that we consider (necessary to understanding 
the bright exciton spectrum) transforms as a 15 dimen- 
sional representation (the adjoint of SU(4)). Such large 
degeneracies of the dark excitons will of course be broken 
if we consider the full forward scattering and/or include 
back scattering terms. But these interactions are rela- 
tively weak in the dark exciton sector and so the breaking 
will be comparatively small. 



Details of Analysis of NRG Data 

In this section we give some of the details surrounding 
the interpretation of the NRG data. In particular we will 
give some details on how finite size corrections determine 
the data analysis. For A^^'^ and A^"''^, finite size correc- 
tions are large because both of these excitations arc near 
two-exciton thresholds. 



But we first turn to the easier of the gaps to deter- 
mine: A-^"'-^, A-^"'^, and A'''""'^. For specificity we focus 
on A^"'^ . Presented in Fig. 4 is A^"'^ as a function of sys- 
tem size, Rnrg J for two different truncations of the total 
number of states (15000 and 75000). (These truncations 
are determined by truncating the system so that the en- 
ergies of the unperturbed states satisfy RnrgE /2tt < 12 
and RnrqE /2Tr < 14). We see that A^"'^ possesses only 
a weak dependence on Rnrg- This, however still leaves 
the question of which value of Rnrg to use for A^"'^. For 
small values of Rnrg, we expect the energies of the full 
theory to reduce to that of the theory's UV fixed point, 
i.e. we expect the energies to behave as 1/Rnrg- We 
can see this behaviour beginning already at Rnrg = Stt. 
For large values of Rnrg we expect instead a truncation 
error (coming from ignoring higher energy states) of the 



form R 



NRG 



This error typically leads to an over- 



estimate of the true value of the gap. The slow increase 
of A-'^"'^ that begins at Rnrg = IOtt is evidence of this. 
Given these behaviors (deviations of A^"'^ from its true 
value), it then makes sense to choose Rnrg where A^"'^ 
is at its minimum. This we do by choosing Rnrg = Qtt. 

The one cautionary note we make is that at this value 
oi Rnrg, A^"'^ may still see finite size corrections. Even 
though Rnrg is sufficiently large so that theory is far 
enough away from its UV fixed point, A-'^"'-'^ can see cor- 
rections related to virtual processes by which the excita- 
tion emits a virtual particle which then "travels around 
the world" and is reabsorbed. Such processes are sup- 
pressed exponentially in Rnrg leading to small correc- 
tions to Ai"'\ i.e. (5Ai"'i - e~"^. The coefficient a is 
either equal to A/vchar where A is one of the other gaps 
in the theory and Vchar is its characteristic velocity (a so- 
called F-term [25|) or it is related to the binding energy 
associated with thinking of A^"'^ as a bound state of two 
other excitations (a so-called /x-term [25|). However for 
A^"'^, a (whether arising from an F-term or a /x-term) is 
sufficiently large so that the correction should be small 
(it is difficult to make an absolute statement about this 
as we cannot easily estimate the coefficient in front of the 
exponential which in principle could be large) . Believing 
it to be small we make no attempt directly to correct our 
results for this error. However dealing with finite size 
errors is necessary in studying the gaps A^^'^ and A-'^"''^. 

We do however do a renormalization group improve- 
ment (RGI) of our numerical data by taking into account 
the RG flow that is occurring at fixed Rnrg ■ We see that 
as we increase the cutoff energy, Ec, the value of A^"'^ 
flows to smaller values. This flow, at lowest order, is 
governed by the equation 
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where /? here is a function of the anomalous dimension 
of Hgap and (5Ai"^i = A^^^'^iE^) - A^"'^{Ec = 00). Inte- 
grating this equation, we can express the value of the gap 
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2/\^dark ^^ 2A^"'^. In such circumstances strong finite 
size corrections are expected. These corrections take the 
form Sfs^'^^'^ ^ cxp{—aRNRG) where again a is related 
to the binding energy of the exciton. If A^^'^ is a bound 
state of the two dark excitons, a is given by 



a^ = (-2A''°'-'=- A^ff'i) 



2vr 



where Vq appears because the velocity of the dark exci- 
tons is very close to vq as they lie almost entirely in the 
spin/flavour sector of the theory. If instead A^^'^ is a 
bound state of two A^"'^, we have instead 



FIG. 5: The value of A^"'^ for Kc+ = 1/5 as function of the 
system size, Rnrg , and the cutoff energy. 



at finite truncation energy to the gap at infinite trunca- 
tion energy: 
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where 7 is some constant. We thus can use the data for 
two different values of Ec to estimate 7 and so extract 
the value of A^"'^ without cutoff. The results of doing 
so arc plotted in Fig. 4 as A^"'^ (RG improved). 

To arrive at the value of A-'^"'-'^ actually reported in Fig. 
2 of the manuscript, we assume the RG improved value 
represents an upper bound to the value of A^"'^ while 
the value of A-'^"'-'^ at EcR/2Tr = 14 represents a lower 
bound. It is then the average of these two values that is 
reported with the error bar being 1/2 the difference. 

While A^"'^ does not suffer from significant finite size 
corrections, A^^^^ docs. We show in Fig. 5 A^^'^ as 
a function of Rnbg- We see that unlike A-^"'-^, A^^-^ 
steadily increases as a function of Rnrg- Wc understand 
this to be a consequence of A^^'^ being a bound state of 
either two dark excitons, A''"'"'^, or two one-photon ex- 
citons, A^"'^, and that furthermore lies just below these 
two exciton thresholds, i.c A^^'^ is only slightly less than 
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where v here is approximately Vc+ but will not be exactly 
that as the exciton A^"'^ is influenced by both the charge 
and spin sectors of the theory. Thus if A^^'^ is close to 
either of the two exciton thresholds (or if v is large which 
it likely is), a will be small and the finite size corrections, 
even though exponential, will decay only slowly as Rnb,G 
grows. 

To estimate A^^'^, we assume the data to have the 
form A29'i(i?jvj^G) = A2s.i(i? = 00) -h /3exp(-ai?ArflG). 
Fitting the data to this form allows us to estimate 
A'^S'^{Rnfig = 00). We take this as an upper bound 
on the value of A'^S'^{Rnj^g = 00) as we expect that 
at larger RnrGj A^^'^ is overestimated in the numerical 
data because of a finite cutoff, Ec- To find a lower bound, 
we examine A^^'-'^ as a function of momentum. Plotted 
in Fig. 5 is the numerical data for A'^^-^{p, Rnbg) where 
the momentum p is related to Rnrg via 2'it/Rnr.g- 
Where A29'i(p, i?^,^^) and A'^aA(p = 0,i?ArfiG) meet 
then serves as a lower bound for A^^'^. (Incidentally 
the manner of this crossing indicates that the dispersion 
of A^S'^ at small momenta is hole-like, i.e. A^S'^(p) ~ 
^23:1 _p2^2me//.) The value for A'^^A that we plot in 
Fig. 5 of the manuscript is an average of these upper 
and lower bounds, while we take as an uncertainty 1/2 
the difference of the two. 



